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We examine the quantum states produced through parametric amplification with internal quantum 
noise. The internal diftusion arises by coupling both modes of light to a reservoir for the duration of 
the interaction time. The Wigner function for the difi^used two-mode squeezed state is calculated. 
The nonlocality, separability, and purity of these quantum states of light are discussed. In addition, 
we conclude by studying the nonlocality of two other continuous variable states: the Werner state 
and the phase-diffused state for two light modes. 
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I. INTRODUCTION 

The two-mode squeezed state is the paradigm of the 
EPR state for light modes. In recent years such states 
have been experimentally realized and thus have been 
used in a number of applications I*]. In particular, these 
states can be used to demonstrate that quantum mechan- 
ics is nonlocal. Using the polarization states of the two 
modes, a violation of Bell's inequality has been experi- 
mentally realized j2j. In addition, an experiment based 
on parity measurements to demonstrate nonlocality has 
been proposed Q- Using parity considerations, a posi- 
tive everywhere Wigner function has been shown to be 
nonlocal. 

In particular, by using a phase space representation, 
nonlocality for the continuous variable squeezed state 
can be analyzed. Because squeezed states can be ex- 
perimentally demonstrated, a lot of attention has been 
given to studying how noise affects these states. Up to 
the present, nonlocality for a pure two- mode squeezed 
state coupled to an external reservoir has been investi- 
gated This noise has been introduced into the state 
by coupling each mode of the squeezed state to inde- 
pendent external reservoirs. Alternatively, this can be 
interpreted as transmitting the two modes through some 
noisy quantum channel. 

Finding an exact analytic expression for the Wigner 
function is often an operose task. Gaussian Wigner func- 
tions are prevalent because they are exact solutions to 
equations which can be easily dealt with. However, in 
most cases only numerical solutions exist. And in some 
cases analytic solutions exist, but only under certain ap- 
proximations. An exact, non-Gaussian solution for the 
Wigner function for a nondegenerate parametric oscil- 
lator has been obtained for the steady-state Al- 
though the state is a two-mode squeezed state with in- 
ternal quantum noise, the steady state solution exhibits 
no nonlocal features and, therefore, is not useful for tests 
of Bell- type inequahties. 

Parametric amplification generates the two-mode 



squeezed state via nonlinear interactions in a crystal. A 
laser beam passes through a crystal for some time and 
the output is two modes of light which are correlated. 
The degree of correlation depends on the interaction time 
as well as the strength of the nonlinearity. Ideally, this 
process generates a pure quantum state which is non- 
separable. We will consider the nonideal process where 
quantum noise is present inside the crystal. 

The purpose of this paper is to study the nonlocal fea- 
tures of the two-mode squeezed state with internal noise. 
We solve the Wigner function exactly to analyze the dy- 
namics of nonlocality and study the many features of this 
quantum state, such as separability and the behavior in 
the steady state. We begin this paper by considering 
two modes which are coupled via a diffusive nonlinear 
crystal. The Hamiltonian has an interaction term plus 
noise terms which couple each mode to a heat bath. A 
linear quantum Fokker-Planck equation is solved for the 
Gaussian Wigner function of the state. An analysis of the 
steady state conditions is given. In Section III conditions 
for purity and separability of the state are provided. A 
discussion of the class of nonseparable mixed states which 
exhibit quantum nonlocality follows. An application to 
Bell's inequality is given in Section IV, where the effect 
of noise on the nonlocality of the state is discussed. With 
Section V we conclude by examining two types of mixed 
entangled states: the continuous variable Werner state 
and a phase-diffused state. 

II. WIGNER FUNCTION. 

The Wigner function for a two-mode squeezed state is 
well-known to be 
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W2mss{ai,a2) = ^exp[-2cosh(2r)(|aip -t- |a2|2) (1) 

+2 sinh(2r)(Q;iQ;2 + a^a2)] 

for two coherent modes of light ai and a2- The amount 
of squeezing of the state is determined by the parameter 
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r which depends on the nonHnearity of the crystal as well 
as the interaction time of the light propagating through 
the crystal. 

In this paper, we derive the Wigner function for a two- 
mode squeezed state which has internal noise. In this case 
a pure squeezed state is not produced. Rather, there is a 
quantum diffusion process present during the generation 
of the squeezed light. This is intrinsic quantum noise 
which is present for the duration of the interaction time. 

The Hamiltonian, in the interaction picture, describing 
the process of parametric amplification in the presence of 
noise is 



H 



ihK(a\a\ - aia2) + T^{a\t i + aiVl) (2) 

j=l,2 



where k describes the nonlinearity of the crystal. The pa- 
rameter r is a reservoir operator which introduces quan- 
tum white noise into the system characterized by the 
mean photon number n. 

The equation of motion for the density operator which 
describes the quantum state of the two light modes is 
given by the master equation 



dp 
'dt 



K{a\a\p — pd\a}2 — aia2/5 + pdia2) 

1=1,2 

'^ni{2a\pdi - aid\p - puidl) 



(3) 



where (3 = {j3i, P2, P2) ^ four- vector and the oper- 
ator D is the displacement operator for the two modes 



D = D1D2 



(5) 



The matrix V is a 4x4 covariance matrix for the two 
modes. 

A double Fourier transform of the characteristic func- 
tion defines the Wigner function W{<y) for the two 
modes. Using standard procedures the master equa- 
tion © can be mapped into the following Fokker-Planck 
equation for the Wigner function 
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in terms of the real position and momentum variables. 
The transform is given by ai = xi + 1x2 and a2 = 
X3 + ix4, so that xi and X3 are position variables and 
X2 and X4 are momentum variables. The drift matrix, 
Aij, and the diffusion matrix, Dy , are constant matri- 
ces, thus defining a quantum Ornstein-Uhlenbeck pro- 
cess. The drift matrix is 
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and the positive-definite diffusion matrix is 



which is obtained by averaging over the reservoir vari- 
ables. The two field modes are coupled to a bath which 
has a mean photon number given by h and we will as- 
sume that the single photon loss rate for each mode is 
equal so that 71 =72. 

The method to convert the above operator equation 
into a c-number equation is straightforward with the use 
of the characteristic function for the Wigner representa- 
tion defined as 



x(/3) = Tr{Dp) = e-3/3^V/3 



(4) 



D = 4(2n+l) 
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The general solution to a linear, multi-dimensional 
Fokker-Planck equation is known and can be solved ex- 
actly by the method of Fourier transform The gen- 
eral Green function solution for a linear Fokker-Planck 
equation is a conditional distribution which is a multi- 
dimensional Gaussian having the form 



VF(x,t|x',0) 



y/detQ{t) 



cxp 



-(x-e^V)TQ-i(i)(x-e^*x') 



(9) 



r 



The matrix Q is a time-dependent matrix which depends 
on the elements of the diffusion matrix and the eigenval- 
ues of the drift matrix. 

From this solution the unconditional distribution is 
found through 



W{x,t)= J iy(x,t|x',0)iy(x',0)dx'. 



(10) 



The initial condition is taken to be the two-mode vacuum 
state given by 



(11) 



in terms of the real variables. After integration over the 
primed variables and transforming back to the complex 
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variables, we have the following form for the Wigner func- 
tion 



III. PURITY AND SEPARABILITY OF 
GAUSSIAN STATES 



1 / f{ai,a2,t) 
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where 



(12) 



}{ai,a2,t) = ci{t)[aia\ + a2a2) (13) 
+ C2{t){aia2 + OL\a2) 



A. Purity 

In terms of the Wigner representation, a general two- 
mode Gaussian state may be written as 



W(a.) 



VdetW i^iW 



e 2- 
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with time-dependent functions 

ci = 2(e-P^ -He-Pi) 

/I — e^P^ 1 - e^P2 

-h (2n-M)(pi +P2) + 

V Pi P2 

C2 = -2{e-'P'' - e-P') (14) 



+ (2n+ l)(pi +P2) 



Pi 



P2 



and 
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+ {2n + 1) 
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2 y V Pi 
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The dimensionless time parameters pi = d + 2r and 
P2 = d — 2r have been used. We call d — "ft the dif- 
fusion parameter and r = nd the squeezing parameter. 
It should be emphasized that ci = ci(t),C2 = C2(i), and 
h — h(t) and, throughout this paper, we occasionally 
omit the function's dependent variables for simplicity. 

The eigenvalues of the drift matrix are doubly degen- 
erate and are 



^1.2 = 
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•^3,4 = - ^(7 
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The condition for a steady state solution to exist is that 
the eigenvalues of the drift matrix all have negative real 
parts. Thus, the parameters pi and p2 alone determine 
whether a steady state solution for the Wigner function 
exists. If ^ > 1 then a steady-state solution exists. The 
steady state Wigner function is a squeezed thermal state 
given by 



W 



exp 



- (|ai| 



\a2\ 



^ (aiaa + a\a*2) 



(V2)^(2n+l)V l-( 



(17) 

This steady state solution is undefined if 2k ~ 7. If 
K = the steady state solution is a thermal state. For the 
case of a pure two- mode squeezed state, 7 = 0, there is 
no steady state solution because as the interaction time 
approaches infinity the Wigner function approaches an 
EPR state with perfect correlations. 



For the state in Eq. (|12|) . W takes the simple form 
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J_ I Ci C2 

2/i I C2 ci 

C2 Ci 



(19) 



The V matrix can be obtained from the W matrix 
through the relation between the Weyl- Wigner character- 
istic function and the Wigner function. They are related 
by a Fourier transform which leads to 



W = EV^E 

where the matrix E=diag[l, -1,1,-1]. 
that 



(20) 



From this we have 
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VdetW 
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so that both matrices differ only by a factor. 
The form of the covariance matrix is 
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(22) 



with the identification + 1/2 = 2hci/{c\ — c^) and 
M = 2hc2/{c{ — c\). N quantifies the correlation {d\di) 
for mode i, while M quantitifics the correlation between 
modes related to squeezing (0,10,2). 

Gaussian operators which are projectors represent pure 
states . The condition for a Gaussian Wigner function 
to represent a pure state may be written concisely as 
vdetW = 4. The Gaussian operator corresponding to 
Eg. lfT^ is a projector provided 



Ci(tf -C2{tf ^l&h{tf 



(23) 



From this condition, it is established that only one pure 
state exists, namely, that for which both r 7^ and 7 = 0. 
This corresponds to a squeezed state which is generated 
in the absence of noise. The presence of the internal 
quantum noise always produces a mixed state. The re- 
sulting mixed states may be separable or nonseparable, 
depending on the parameters of the model. 
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B. Separability 

A separability criterion for two-mode Gaussian states 
has been established which relies on the partial trans- 
position map acting on the two-party state @. It has 
been shown that this criterion is equivalent to determin- 
ing that the quantum state is P-representable. A two- 
mode Gaussian state is P-representable, and hence, sep- 
arable, if and only if 

V-il>0 (24) 

where V is the covariance matrix for the characteristic 
function found in Eq.Q) which takes the special form 
of Ea. (|22|l . The separability condition is also equivalent 
to the condition that the matrix elements of (|22|l satisfy 
iV > M. 

The state in Eq. H12|l is separable if and only if the 
eigenvalues of Ea. (|24|l are nonnegative. There are four 
eigenvalues which are doubly degenerate. They are 
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FIG. 1: A parameteric plot shows the states which are sepa- 
rable. States which lie on or above the solid line are separable 
states. The dash-dotted line is for d = 5 while the dashed line 
is for d — 2.5. The two dots correspond to n=0. The mean 
photon number increases from to 10 along both lines. The 
squeezing parameter is fixed at r — 1.5. 



IV. NONSEPARABLE MIXED STATES 
EXHIBITING QUANTUM NONLOCALITY 



(1 - e-P''){dn + r) 
ei,2 = 

P2 

(l-e-PO(dn-r) 

63,4 = • (25) 

Pi 

The first eigenvalue has the property that the sign(ei,2) 
is positive for all parameter space. It is not important for 
establishing the nonseparability of the state. The second 
eigenvalue determines the nonseparability of the state. 
For certain regions of the parameter space it is negative, 
while in other regions it is positive. The negativity of 63^4 
may be summarized as sign(e3,4)=sign(c?n — r). Hence, 
the state is nonseparable if and only if 

r > dn. (26) 

This is the requirement that the squeezing parameter be 
greater than the product of the single-photon loss rate de- 
scribed by the internal diffusion parameter and the mean 
photon number of the reservoir. Fig. 1 shows a paramet- 
ric plot for varying mean photon number n. States which 
lie on or above the solid line are separable states. The 
dots show states which have n = for different values 
of the diffusion parameter d. We immediately see that 
if n = then 62 < for all parameter space so that the 
system is never separable. Therefore, one may conclude 
that this system is entangled for all of parameter space 
when coupled to a zero-tempcrature heat bath. This is 
interesting because the diffusion parameter may be made 
as large as desired and still the state is nonseparable, 
even as the diffusion approaches infinity. The nonsepa- 
rable states may be further classified into two sets: those 
which exhibit quantum nonlocality and those which do 
not. 



If the two-mode squeezed state is generated with inter- 
nal noise present then the resulting state will be a mixed 
state. Although all pure entangled states violate some 
Bell inequality, it is not clear in general which mixed 
entangled states will do so. If no noise is present then 
there always exists, for a fixed r, some range of values 
for J which will correspond to a nonlocal state. The 
presence of the diffusion parameter will destroy the non- 
local features until some critical value for d is reached 
such that the state goes from nonlocal to classically cor- 
related. Thus, in addition to the pure states {d — 0) 
which are nonlocal, there is a set of mixed states {d 7^ 0) 
which is also nonlocal. 

To study the nonlocal properties of the state l|12|) 
we use the formalism of parity operator correlations 
developed in The Wigner function is related to 

the expectation value of a displaced parity operator in 
phase space given by n(ai, 02) = (-D(ai)(— I)"^i)^(ai)(8) 
^(a2)(— 1)"^-D^(q!2))- This expectation value defines a 
correlation function 

n(ai,a2) - (I) V(ai,a2) (27) 

for parity measurements of mode I and mode 2, which is 
proportional to the two-mode Wigner function. In anal- 
ogy with the case of dichotomic spin states, a Bell combi- 
nation may be written, using correlations between parity 
measurements, as 

B = n(o,o)+n(%/j,o)-i-n(o,-\/j)-n(yj,-Vj). (28) 

For a two-mode squeezed state with internal noise the 
combination B will be a function of the interaction time 
as well as the squeezing and noise parameters. The Bell 
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combination is 



B{J,t) 



1 

W) 



1 + 2 exp -J 



2h{t) 



(29) 



exp ~ J 



h{t) 



where the functions ci, C2, and h are given by Eas. (|14ll5|l 
which are imphcitly functions of time, and J is the square 
of the coherent amphtude of the Ught mode. 

A search algorithm finds a local maximum for B and 
returns the values of the free parameters. The function B 
has a maximum of 2.19 with values J = .01, r = 1.5, d = 
0, and n — Q. As expected, the Bell combination is most 
nonlocal when there is no noise present and the squeezed 
state is a pure two-mode squeezed state. In this case, one 
finds a stronger violation for large r and small intensity 
J. To see how the internal noise affects the nonlocality of 
the state we set r = 1.5 and let J and d vary in Fig. 2. It is 
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FIG. 3: The Bell combination is plotted as a function of the 
diffusion parameter, d, with the squeezing parameter fixed at 
r = 1.5, n = 0, and J = .01 



state, and a convex combination of a pure entangled state 
with a classically correlated state. The first case is the 
analogue of the Werner state for continuous variables. 
Werner |0 investigated the state 



1 - n » 

p = p|*)(*| + — 2^/a«)/b 



(30) 



for < p < 1 which is a convex combination of a maxi- 
mally entangled state l^*) = X^iLi \''')A®\i) b for finite 

dimension d^ with a maximally mixed state obtained by 
a partial trace over the maximally entangled state. It 
was shown that such a state is nonseparable if and only 
if 



. 04'- 



1 



(31) 



FIG. 2: The Bell combination is plotted as a function of the 
coherent amphtude, J, for the two field modes and the dif- 
fusion parameter, d, with the squeezing parameter fixed at 
r = 1.5 and n = 0. 

clear that the diffusion parameter destroys the nonlocal 
features of the state. Fig. 3 shows a cross-section of 
the three dimensional graph of Fig. 2 for J = .01. The 
Bell combination is greater than 2 for small values of the 
diffusion parameter, d. On a large scale, we find that 
as d increases the Bell combination decreases, reaches a 
minimum, and then increases to approach 2 from below. 



V. CONTINUOUS VARIABLE WERNER 
STATES 

Using the two-mode squeezed state of Eq.Q, a class 
of nonseparable, non-Gaussian, mixed states can be ex- 
plored. We will consider two cases: a convex combina- 
tion of a pure entangled state with an uncorrelated mixed 



A. Nonlocality of the continuous variable Werner 
state 

A natural extension of a Werner state for infinite di- 
mensions is a convex combination of a maximally entan- 
gled state with the maximally mixed state 

PW = pptls + (1 - p)Pt <» Pt (32) 

for two modes of light, A and B. This state is analogous 
to the Werner state when the squeezing parameter r is 
infinite. The continuous variable phase-space represen- 
tation of this state is given by 

Ww = pW2mss{ai,a2) + (1 - p)WT{ai)WT{a2) (33) 

where Wt{o:) is the marginal Wigner function obtained 
by the integration 

Wria) = / W2mss{ai,a2)d^a2 = 7^"! (34) 
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over the variables of one mode. When r is infinite, the 
Wigner function W2mss is the EPR state and the function 
Wt{o:) has infinite variance representing a maximally 
mixed state of knowledge. Note that Ea. (|33|l is not a 
Gaussian operator, but rather, it is a linear superposi- 
tion of Gaussian operators. The separability conditions 
for a linear superposition of two Gaussian operators have 
not yet been established. However, it has been shown 
that the state pw is nonseparable if p > ^ . This is 
the limit of the finite dimensional case, Ea. (|31|l . as the 
dimension of the Hilbert space, d, approaches infinity. 

A test of nonlocality for the state Ww in phase space 
can be performed using the same Bell combination given 
by Ea. H28|l . It is found that the state remains nonlocal for 
a region of mixtures lying between p = .9 and 1. Adding 
a mixed state to the entangled state rapidly degrades the 
nonlocal features. 




given by 

W= 4e~'™'''^''^*'"'''+'"'''^/o(4|ai||a2|sinh(2r)) 

(37) 

where Iq is the modified Bessel function of the first kind. 
In contrast to the Werner state, this state does not fac- 
torize into a product of functions containing each mode. 
The state ly is a classically correlated state. The state 
of Eq.JSnj now being considered is a mixture of a state 
containing quantum correlations and a state containing 
classical correlations. The nonlocality of this state is ex- 
amined by considering an expansion for small intensity, 
J, of the Bell combination given by Ea. (|28l) . This expan- 



(38) 



B = 2 + 4Jpsinh(2s) + 0{J^ 



reveals that the state (|35ll is nonlocal for all p > 0. This 
is seen in Fig. (5) which shows that only the case p = 
washes out the nonlocal features completely. Thus, 
we have presented a classically correlated state which, 
when mixed with even the smallest amount of a quantum 
correlated state, shows nonlocal features. 
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FIG. 4: The Bell combination is plotted for the continuous 
variable Werner state. The squeezing parameter is fixed at 
r — 1.5. The curves are shown for the following values of 
mixture, p, from highest maximum value to lowest maximum 
value: p — 1 (top-most curve), p = .95, p = .9,p = .5, and 
p = (bottom-most curve). 



B. Nonlocality of the phase-diffused state 

In addition to the Werner state of Eq. 1)32(1 for contin- 
uous variables, a different convex combination has some 
interesting properties. Let us consider the state 



W ^ pW2mss + (l - P)W 



(35) 



which is a convex combination of a two-mode squeezed 
state with a phase-averaged state. Rather than integrat- 
ing over one mode of light, an integration over the phases 
of both modes is performed. The state 



W = 



W2n 



2tt 2tt 



(36) 



describes phase diffusion of two light modes such that 
each mode has a completely random phase. The inte- 
gral of Eq. H36|l leads to a phase-averaged Wigner function 
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FIG. 5: The Bell combination is plotted for the Werner-type 
state pW2mss + (1 — p)W. The squeezing parameter is fixed 
at r = 1.5. The top-most curve shows the pure two-mode 
squeezed state obtained for p = 1. The remaining curves, 
from the greatest maximum value to the lowest maximum 
value, show p — .5,p — .2, and p = (no violation). 



VI. CONCLUSION 

In this paper we examined the two-mode squeezed 
state with internal quantum noise. The Fokker-Planck 
equation was solved to obtain an exact solution for the 
Wigner function. The steady-state and nonstcady-state 
regimes were explored. We found that nonlocal features 
of this state exist in the regime where no steady-state 
solution exists. 

The parameter space contains separable and nonsepa- 
rable states, both mixed and pure. Only one pure state 
exists which is nonseparable. Nearly all states in the 
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parameter space are mixed states. Using the Gaussian 
Wigner function for the state, the separability criterion 
was determined to identify those mixed states which 
are nonseparable. There exists a set of nonseparable 
mixed states for which the quantum diffusion can be large 
< d < oo. 

Using the Wigner function, a test of nonlocality in 
phase space determined the effect of internal quantum 
noise on Bell's inequality. As expected, the noise param- 
eters reduce the nonlocality of the state. But there is 
still a region of mixed entangled states which are non- 
local. To further study nonseparable mixed states ex- 
hibiting nonlocality, wc investigated Werner-type states 
for continuous variables. Wc found that mixing the two- 
mode squeezed state with a product of two thermal states 
destroys the nonlocal features state. Additionally, we 



found that mixing the pure two-mode squeezed state with 
a phase-diffused, classically correlated state did not de- 
stroy the nonlocal features, regardless of the amount of 
mixing. With this, we have provided an example of a 
classically correlated state which, when mixed with even 
the smallest amount of pure entangled state, is nonlocal. 
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